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SUMMARY 


This  report  covers  the  first  eighteen  months 

of  our  investigations.  Following  a  brief  • 

introduction  the  report  discusses  in  some 
detail  the  mathematical  model  of  the  system. 

the  methods  used  in  solving  the  corresponding  • 

differential  equations  and  the  various  ways 
of  presentation  of  the  resilts.  The  effect 

of  an  external  load  is  then  discussed  in  some  • 

detail.  The  report  ends  with  suggestions  for 
future  work  based  on  the  results  obtained 

so  far.  # 


-  5  - 


1.  Introduct ion 

The  main  purpose  of  this  work  is  to  investigate 

the  control  of  chaos  in  a  plasma.fi 1 led  diode. 

A  mathematical  model  of  such  a  system  has  been 

1  ° 

considered  before  '  ,  its  main  feature  being  the 
replacement  of  a  smooth,  interelectrode  electron  cloud 
by  a  large  number  of  infinitely  thin  electron  sheets. 

In  the  absence  of  ions  similar,  continuous  systems 

Tt 

have  been  considered  in  great  detail  '  ,  including  the 

effects  associated  with  the  granular  nature  of  the 

electron  stream  .  In  the  presence  of  positive  ions 

the  system  was  first  considered  by  J.R. Fierce  \ 

6  7 

and  then  in  more  detail  by  others  ’ f  .  When  the 

electron  cloud  is  discretized,  the  ions  can  be  treated 

as  a  smooth  positive  background  or  they  can  be 

discretized  ns  well.  Drawing  on  the  experience  quoted 
!* 

above  we  are  suggesting  that  the  latter  is  the  more 
appropriate  approach  to  adopt  for  further  investigations. 


...  • 

. —  -wgy*  •  ^ 


o 


Mathematical  model 


Let  us  first  of  all  consider  a  system  which  is  devoid 
of  positive  ions  (pure  vacuum).  The  equation  of  motion 
for  a  single  sheet  of  electrons  is  then  given  by 


(2.1) 


where  x.  is  the  position  of  the  sheet,  p  and  p 

1  *  is  '  ms 

are  the  respective  charge  and  mass  surface  densities 

e/m^)  and  ®  is  the  average  field  acting 
on  the  sheet. 


t(*T 


♦ 


f? 

TI 


) 


( 2 .  :> ) 


Here  and  ETT  are  respectively  the  fields  in  the 

grid/electron-sheet  and  electron-sheet/anode  spaces, 
the  general  geometry  being  that  shown  in  Tig.l.  Since 
the  field  must  change  from  to  ns  ve  Pass 

through  the  sheet,  the  act  lal  fieli  acting  on  the  sheet 
is  the  average  field  v  .  Substititing  (2.2)  in  (2.1) 
we  find  th  t  the  single  sheet  Is  decelerated  when  x^  <  d/2 
and  accelerated  when  x^  >  d/2. 


«•  t» 


g 


7  - 


Tn  the  cnse  of  several  sheets  the  nver.^e  field 
acting  on  sheet  j  is  qiven  by 


s. 

J 


a  sheets  to  the  left  of  j  )  ♦  S  ^  ♦ 

KT(all  sheets  to  the  right  of  j) 

-  ^TI(i>  j)  ♦  B  *  Bx(i<  j)  .  > 

•  ~t[  &"*•  *  -  i)  *  £  *H 

■  -ft  [  &  -fr  -  *  *  i] 


(2.?) 


where  the  sheets  are  counter!  from  ri  jht  to  left,  since 
sheet  l  would  he  injected  first  and  would  have  travelled 
to  the  right  by  the  time  sheet  2  has  been  injected,  and 
so  on.  Here  M  represents  the  total  njmber  of  sheets 
momentarily  in  the  interelectrode  space. 

In  ilace  of  (2.1)  we  now  have  the  following  set 
of  M  coupled  equations  of  motion,  each  representing 
one  of  the  M  electron  sheets  present  in  the  interelectro  Ip 
space : 


•  • 

1  x, 
ms  1 

■  e.\ 

« 

4 

• 

# 

» 

t 

>  T  . 

ms  j 

"  Ps?j 

• 

0 

• 

» 

# 

4 

’msXM 

=  Ps?M 

(2.4) 


It  is  now  convenient  to  introduce  the  following 
reduced  variables  *  : 


x/d 


t/t 


(2.5) 


where  d  is  the  separation  of  the  electrodes  and 


t 

o 

is  the  transit  time  in  the  nbsenca  ©f  fields  of  an  electron 

Injected  with  initial  velocity  x  =  v  ,  i .  t  -  d/v  . 

o  o  o  o 

Thus  we  now  have : 


dx  1  t»  dx  x  d*  x  1  t*  d*  x 

dt  '  *  d  dt  vD  *  dt  '*  "  d  dt* 


In  order  to  introduce  the  concept  of  convection  current 

J  we  write 
c 


D  =  -ffjU-T’T..  »  .  _ik»_d_ 

rs  n5  v0  x 


(2.7) 


where  is  again  the  ehar/?e  surf ace- densi ty  of  a  single 

sheet,  p  is  the  surface  density  of  X  sheets  and 

S  ,  l  O  t 

po  is  the  smoothed-out  volume  density  of  X  sheets. 

\  typical  equation  (2.4)  can  now  be  written,  substituting 
from  (2.7),  (2.6)  and  (2.?): 


di 

dt 


x1; 

757" 


P>  _  J**  t® 

pn*  ‘tV.V 


(2.8) 


Xote  that  M  may  be  different  f rom  X  ,  since  X  is 
a  constant  and  M  may  vary,  being  equal  to  the  number 
of  sheets  which  are  present  in  the  interelectrode  space 
at  any  given  time. 

In  the  final  step  of  our  derivation  of  the  eq  tations 
of  motion  we  choose  to  express  the  convection  current  -T 

co 

as  a  fraction  of  the  convection  current  flowing  in  a  similar 


sp^ce-chnrge-1 imi ted  diode  T  given  by: 

c  •  s pi 


t  ,  =Jh.  _£a>  »  ilst.  -£a»  -\- 

Jc,scl  9  p*  rt*  9  pj 


(2.9) 


Substituting  (2.9)  in  (2.8)  we  finally  obtain 


-fr  [ I,  s'i  -  j  *  }. 


(2.10) 


t  •  J  /J 

CO  c  ,  scl 


(2.11  ) 


ns  a  typical  equation  of  notion  for  the  jth  electron 
sheet  1  . 

The  M  equations  (2.10)  fully  describe  the  motion 

of  M  electron  sheets,  each  sheet  injected  with  the  same 

initial  velocity  dx  '  . /dt 1  *  x  /v  ■  1  into  the  interelectrode 

J  o'  o 

space.  Since  the  diode  is  assumed  to  be  short-circuited, 

the  only  fields  present  are  those  due  to  the  space  charge 

of  the  sheets  themselves.  Thus  (2.10)  represent 

a  computer  model  of  the  conditions  inside  a  plane  diode, 

where  the  cathode  has  been  replaced  by  a  grid,  so  that 

the  electrons  can  be  injected  with  a  constant  velocity 

3 

v  which  is  different  from  zero  . 
o 

So  far  we  have  only  considered  a  short-circuited 
diode.  If  a  load  is  connected  across  the  electrodes, 
we  require  an  additional  equation  which  relates  the 
conditions  inside  and  outside  of  the  diode.  Voting  the 
fact  that  the  total  current  in  Fig.  1  must  be  continuous 


we  can  write: 


m i  i® 


-10  - 


T  ♦  T 

i  nd  cap 


d  v «  *  l 


(2.12) 


where  V  is  the  potential  difference  across  the  load 
and  a  suitable  expression  for  the  induced  current  I  ^ 
has  been  derived  elsewhere  *  .  In  the  case  of  a  load 
consisting  of  a  resistor  R  in  series  with  an  inductor  L 
we  have 


V  -  RT .  _  +  LI„  _ 

tot  tot 


(2.1?) 


so  that,  substit  iting  in  (2.12)  we  finally  obtain 


■  — C  X-  x .  —  C  RI 


Vi-,  i 


V  tot 


-  C  LI 


v  tot 


(2.l4) 


where 


(2.15) 


expresses  the  plane  (geometry  of  the  system.  One  can  also 
look  at  it  as  an  irreducible  capacitance  of  the  system. 

In  the  presence  of  the  load  the  equations  of  -notion 
(2.4)  must  be  adjusted  by  writing: 


(g.  ♦  e,  „>  -  -e*.  it 

p«i*  j  load  j 


-  j  +  *1  -  t  ♦  Lf  )  (12.16) 

to  Lv»*  d  J  a.  !  tot  tot 


Multiplying  both  sides  of  ( 2 . 1  ’■+ )  by  tQ/psA  and  (2.15) 

o  . 

by  t'  /d  we  obtain  corresponding  expressions  in  terms 
of  the  reduced  variables  (2.5)  -  (2.?): 

dt^r  -  f?T  |,x,i  -  1  *  -  V  -  tl  dt Tl 


(2.?)  : 

T  i »  r  ^ 1 

"  TRX  TL  dt'J 

(2.17) 

dal 

dt  ' 

(2.13) 

Here  we  h;  ve  introduced  a  new  reduced  variable  given  by 


.A  tot 


(2.19) 


The  constants  T_  and  T,  are  defined  by  writing: 
K  L- 


-fit-. -SsL  .  -jgt,. 

d  9N 


(2.20) 


-£*-.-£-ri.  .  -  -  -fit J&A.JL  T 

I °«*  *  tot  te  N  ^  d  d  t4 


P-e  <r 


2t 

C*R 

T  1  - 

2  L 

r  1 

9V 

to 

1  “ 

9\ 

tr  ] 

l 

dT  1 
dt  ' 

- 

P-* 

-r- 

t* 

, -is., 
d 

2  k 

CwL 

dl ' 

2c 

.  T 

dl 1 

9V 

dt  * 

9N 

*  i 

'L 

dt  ' 

(2.21) 


L  dll 

If  dt* 


(2. ?2  ) 


Since  <  0  in  (2.19).  it  should  be  noted  that  I* 

has  an  opposite  sign  to  I  In  the  case  of  Fig.  I  we 

have  I.  .  <  0  and  I*  >  0. 
tot 

One  can  see  from  (2.17),  (2,13)  that  in  the  presence 


of  a  RtL  load  M+l  coupled  second  order  differential 
equations  are  required  for  the  description  of  the  system 


12 


In  some  coses  it  is  preferable  to  use  the  field 

i?  i  d  ,,  ed  _  p*  d 

h  *  »  —gp  ■  e  *  ■  ■  ■  j  w  u  „  r  p 

2%  load  ^  v©  load  p^s  v*  load 

•  •■•*...>  -  «■  * 

*£#^)-^«v  *TtiK 


T5  *  ♦  E « 

R  L 


(2.23) 


where,  from  (2.23) 


_L_  dKW 

T0  dt  ' 


(2.2U) 


Substituting  for  I*  in  (2.17),  (2.13)  we  now  obtain 


it7*  9V  [  &,  i 


i  .  i  »i  T>,  dR,[ 

J  +  tj  ~  E  R  "  TR  dt  ' 


(2.25) 


Jk  Ml  _1_  pi  ,  2C  ^  dx-J 

T*  dt,x  "  dt  '  "  Tr  R  9M  ;-i  dt  1 


(2.26) 


The  evolution  of  the  system  is  again  described  by  M+l 
coupled  second-order  differential  equations. 

Another  load  to  be  considered  is  in  the  form  of  a  resistor  R 
and  a  capacitor  C  in  parallel.  In  place  of  (2.12),  (2.13 
we  now  have 


*  tot  "  ^ind  *  ^cap 


-  -ifit-  f  x  -  C  V 
d  xi  "vV 


-  TR  +  xc 


+  cv 


(2.27) 


W  19 


J 


8 


-  1>  - 


where  I_  and  I  respectively  represent  the  currents 
H  C 

flowing  in  the  two  branches  of  the  load.  Since  the  field 


across  the  load  is  given  by  the  usual  expression 
-  V/d,  we  obtain  from  (2.27) 


"load 


B 


load 


-?-*¥■ 


•-e-% 


-  -£l  t 

U  'R 


CvR 


dx‘J 
dt  1 


(tr  *  V  -Tt*4  <2-23> 


where  we  have  «sed  (2.5),  (2.15)  and  (2.20)  _  (2.22)  foi 

defining  x 1 ,  t1,  T^  ;  in  addition  we  have  introduced 
a  new  constant  given  by 

Tc  -  RCAo 

From  (2.2?)  we  can  now  write  for  the  reduced  field 


(2.29) 


_ 

pm»  V 


p£d*>  Ri'ldt'  ‘k+~C  dt 1 


or 


E  1  ■ 


2  L 
9N 


N 

T  X 

r  r.i 


dxt1 
dt  1 


-  (TS  *  V 


dE  1 
dt 1 


(2.30) 


This  equation  should  be  compared  to  (2.13)  which  is  valid 
for  an  R  and  an  L  in  series.  The  corresponding  equation  of 
motion  of  the  electron  sheets  is  now  given  by 


M 

v  2L 

• 

X  .  t 

1 

K(Cy 

+  C) 

• 

-V 

d 

M 

K~ 

AsL  + 

—  ( c 

V  V 

+  C) 

dEiowi 

• 

»  - 1 

dt '  " 

dt  ' 

ft!  ftft 


d*  X!1 

dt™ 


*i  -  J  *  *1  -  E' 


(2.31) 


This  equation  is  similrr  to  (2.25)  except  for  S'  in 

place  of  the  two  terms  in  S3 . 

K 

Sometimes  it  is  more  convenient  ro  express  everything 
in  term-,  of  current.  7e  then  have 


'load 


-  RVd 


(2.32) 


-  E 

p^d  load 


>jt 

* 


R  ^7 


fcT  Xt,  d  JR  jC&i  Tr 


— -  T  T  • 

9N  R  R 


(2.33) 


Thus  in  terms  of  current  (2.31)  and  (2.30)  acquire  the 
following  algebraic  form: 


d*  x  ’• 


- j  *  *1  -  vr! 


<tr  *  Tc>-^-  *  H 


(2.34) 


(2.35) 


It  should  be  noted  that  in  the  load  equations  the  capacitance 
C  is  merely  added  to  Cy ,  the  two  being  in  parallel 
(  Tr  +  Tc  -  R(Cy  ♦  C)/T.  ). 


V!  ftp 


Let  us  now  consider  briefly  the  initial  conditions. 

Tt  is  assumed  that  initially  the  diode  is  free  of  charge 
and  that  no  external  batteries  are  present,  the  energy 
being  delivered  to  the  system  by  the  injection  current  t 
In  those  circumstances  there  are  no  charges  on  the 
electrodes  until  the  first  electron  sheet  is  injected 
into  the  interelectrode  sp  ce  through  the  left-hand 
electrode,  i.e.  the  grid.  Consequently  in  the  case  of 
(2.17),  (2.13)  or  (2.25),  (2.26)  the  initial  boundary 
conditions  are  respectively  either 

(2.36) 

(2.3 7) 

In  the  case  of  (2.30),  (2.31)  or  (2.34),  (2.35)  we  again 
have  either 


d  T  1 

I  1  -  o,  .0  at  t '  »  0 

or 

E '  -  0 ,  -  0  at  t  '  -  0 

H  clt  * 


r\V  * 

E '  -  0,  ~r  -  0  at  t  '  -  0 

at 

or 

TR  "  °*  dt^  *  0  at  *  0 


(2.33) 

(2.39) 


Finally  let  us  consider  an  approximate  expression 
for  the  potential  distribution  between  the  electrodes. 
Since  the  interelectrode  space-charge  in  our  model  is 
represented  by  a  succession  of  electron  charge  sheets, 
the  potential  function  is  approximated  by  a  series  of 
straight  line  segments  connecting  the  potentials  of 
neighbouring  sheets. 


.  16  . 


Expressions  for  dj  and  are  obtained  by 

integrating  the  corresponding  fields  and  -'jj 

<1  loted  in  (2.2);  since  in  our  model  the  fields  are 
constant,  the  corresponding  potential  functions  must 
by  straight  lines  sloping  downwatds  between  the  grid 
and  point  i  in  Fig.  1,  and  upwards  between  i  and 
the  anode.  For  example: 


*z  -  -  w 


ijs-  ,  <  o 

e*  d 


(2.1*0) 


eti  rtx  *  -  ~t  lx~d)  *  ° 


(2.41 ) 


since  j>s  <  0.  In  the  presence  of  many  sheets  we  find 


that  for  sheet  j  we  have  : 


- .  ?  ^TT1  ♦  z.  *L  + 


TTi  Cl,  T  i 


tSt  i h(d  -  V  *  \j(d  -  hO 
\  hUd  -  I,  V  *  dXxi]  *  +- 


(2.42) 


In  the  presence  of  a  load  (2.29)  has  an  additional 

terms  given  by  V.  -  -  E.  .x, 

j  load  j  .  Dividing  both  sides  by 


-If-? 


(2  .43) 


which  is  another  form  of  (2.9)  and  then  using  (2.7) 


17  - 


we  finally  obtain 


An  approximation  to  the  potential  function  fl 1  is 
obtained  by  connecting  all  the  points  .  with 
straight  line  segments. 

It  should  be  noted  that  (2.17),  ( 2 . 1 S )  and  also 
(2*25),  (2.26)  are  linear;  however  the  boundary 
conditions  of  the  system  are  expressed  in  terms  of 
(2.M)  which  is  nonlinear,  so  that  the  inherent 
nonlinearity  o f  the  system  is  suitably  expressed  in 
this  fashion. 


»«  ttt 
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2* _ Computational  details 

Two  different  methods  have  been  used  for  the 
solution  of  equations  (2.17),  (2.13)  and  (2.25), 

<2.26)  or  <2. ?0),  (2.31)  rnd  (2.3*4).  (2.35). 

The  first  method  is  based  on  the  use  of  backward 
and  forward  differences.  \s  is  customary  in  such  cases 
each  second  order  differential  equation  is  first  converted 
into  two  first  order  differential  equations  and  the 
corresponding  derivatives  are  then  expressed  in  terms 
of  the  finite  differences,  using  second  order  approximation. 
Thus,  for  example,  we  obtain  for  the  equations  of  motion 
of  the  j-th  electron  sheet: 

Xk+1  *  |  2  Xk  ‘  *  Xk-1 1  A  t  +  Xk 

and 

xk  +  l  *  lxk  +  (  f  xk  -  S  xk-l)AtlAt  *  xk 

This  method  has  been  initially  developed  by  Birdsall 
and  Bridges  ^  for  the  solution  of  equations  describing 
the  conditions  in  a  short-circ  lited  diode.  We  have 
extended  it  to  cover  the  case  of  a  diode  with  a  partly 
reactive  load  connected  across  its  electrodes. 

The  second  method  of  solving  the  M  +  1  coupled 
differential  equations  is  based  on  a  fourth-order 
Runge-Kutta  subrojitine  of  proven  reliability.  This 
method  is  particularly  convenient  when  the  number 

of  equations  is  large,  say  greater  than  20;  it  is 
also  very  fast. 


(3.1) 

(3.2) 
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The  following  methods  have  been  used  for  testing 
the  accuracy  of  our  results.  In  the  absence  of  a  load, 
l.e.  when  the  diode  was  short-circuited,  ve  would 
compare  our  numerical  results  with  those  obtained 
analytically .  In  most  cases  the  agreement  was  excellent 
(of  the  order  of  0.01^)  provided  the  injection  c  irrent 
U  <  8;  at  t  3  a  virtual- cathode  appears  and  the 
sheet  model  develops  oscillations. 

For  r>  irely  resistive  loads  we  were  able  to  use 
Birdsall's  and  Runge-Kutta  subroutines  and  then 
compare  the  results;  as  a  rule  we  would  reject  any 
results  that  would  differ  by  more  than  0.l£.  In  the 
case  of  more  complicated  loads  containing  reactive 
components  we  could  only  use  the  Runge-Kutta 
subroutine;  we  would  then  test  the  numerical  accuracy 
of  the  results  by  systems tical ly  reducing  the  time 
integration  interval  until  we  were  satisfied  that  the 
accuracy  was  better  than  0.1^.  This  simple  proced  ire 
can  be  q  lite  effective  in  practice,  provided  it  is  used 
sensibly. 


fit  ftdfe 
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4.  Presentat ion  and  discission  of  the  results 

There  ore  several  ways  of  iresentin^  the  results 
of  our  calculations  and  we  will  discuss  them  in  turn. 

(A)  Sheet  position 

The  simplest  way  of  presenting  our  results  is 

to  plot  the  position  of  all  the  sheets  of  charge ,  as 

shown  in  Pigs.  2  a _rl  for  four  different  values  of  the 

injection  current  density  w  -  J  /j  The  position 

co  s  c  1 

of  the  sheets  is  shown  along  the  horizontal  x'-axis, 
the  vertical  axis  being  used  for  indicating  their 
potential  relative  to  the  entrance  grid  situated  at  x'  »  0. 
Since  in  this  case  the  diode  is  short-ci re  li ted ,  the 
reduced  potential  of  both  electrodes  is  the  same  and 
equal  to  1  -  1.  The  whole  picture  has  been  obtained  for 
a  single  instant  of  time,  well  after  the  transient  due 
to  a  gradual  injection  process  has  died  out.  The  values 
of  the  current  density  c  have  been  chosen  so  as  to  cover 
the  three  areas  of  interest  shown  in  Pig.  3,  i.e.  stable, 
hysteresis  and  oscillatory  (unstable)  regions.  It  is  to  be 
noted  that  in  the  last  two  cases.  Pigs.  2  c  %  d,  there  exists 
a  marked  tendency  for  the  sheets  to  cluster  near  the  > 
injection  electrode,  their  potential  becoming  negative 
in  many  cases. 

(B)  Interelectrode  potential  profile 

Another  simple  presentation  of  the  results  is  that 
shown  in  Pigs.  4  a,b  ,  where  we  have  plotted  the  reduced 
potential  of  each  sheet  against  its  position  along 

the  x*  axis,  all  values  being  obtained  for  a  single  instant 
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of  tin*  and  after  the  transient  has  died  out.  Thus 
the  B-type  presentation  is  very  similar  to  the  A-type 
presentation,  although  now  the  potential  curve  is  much 
easier  to  visualize.  Ye  have  also  used  this  example 
to  show  the  difference  Ln  the  potential  distribution 
between  a  short-circuited  diode.  Fig.  4  a  and  that  with 
a  simple  resistive  load  defined  by  Tn  »  C  R/ 1  »  0.1 

connected  across  its  electrodes.  Fig.  4  b.  Ve  can  see 
that  in  the  latter  case  the  anode  is  at  alower  potential 
than  the  grid,  the  potential  difference  V  »  RI  being 

v  O  v 

generated  by  the  current  flowing  through  the  resistive 
load  R  ,  as  shown  in  Fig.l.  Roth  results  have  been  obtained 
for  the  same  current  density  l  ■  6. 

The  main  weakness  of  the  simple  presentations  of  the  type 
A  and  B  is  that  they  are  valid  for  a  single  instant 
in  time  only,  i.e.  they  do  not  tell  us  anything  about  the 
development  of  the  system  with  time.  In  order  to  overcome 
this  ve  mist  use  other  methods  of  presentation. 

(C)  Tra,jector ies 

Since  the  system  is  IT),  one  can  gain  information  about 

its  evolution  by  plotting  x*  »  x/d  against  t'  »  t/t  . 

o 

Figs.  5  a  and  b  show  the  corresponding  trajectories 
for  V  -  4.0  (short-circuited  diode)  and  i/  *  9.0  (diode 
with  a  R,L  load).  In  the  first  case,  Fig.  5  a  corresponds 
to  a  stable  potential  distribution,  the  charge  sheets  mov€ 
smoothly  between  the  electrodes.  The  second  case,  Tig.  5  b, 
corresponds  to  the  other  extreme  where,  due  to  the  large 
values  of  c  and  we  have  quite  violent  changes  in  the 

sheet  trajectories,  leading  to  large  amplitude  oscillations. 
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Sine*  thf5  position  of  each  charge  sheet  is  noted  at  fixed 
time  intervals  and  then  indicated  by  a  dot,  the  density 
of  the  dots  is  inversely  proportional  to  the  velocity 
of  travel  of  the  sheets. 

(0)  Depth  of  potential  minimum  versus  time 

Strictly  speaking  there  is  no  potential  minimum 
in  our  case;  all  we  have  at  our  disposal  is  a  charge  sheet 
which  is  at  a  potential  that  is  lower  than  the  potential 
of  all  the' other  sheets.  Figs.  6  a  and  b  show  the  value 
of  such  a  potential  as  a  function  of  time  respectively 
for  two  different  values  of  the  in  jection  eurrent  density  i. 
We  find  that  for  i  •  4.0,  Fig.  6  a,  the  lowest  potential, 
after  a  brief  transient,  settles  down  to  a  constant  value. 

In  the  case  of  t-  »  9«0.  Fig.  6  b,  the  lowest  potential 
oscillates  in  a  somewhat  uneven  or  'chaotic'  fashion. 

These  oscillations  have  a  clear  character  of  'relaxation1 
oscillations,  their  upward  movement  being  much  faster 
than  their  downward  movement,  as  can  be  seen  by  the 
density  of  dots.  Closer  investigation  of  drawings  such  as 
Figs.  2  c  A  d  indicate  that  the  oscillations  occur  because 
the  charge  sheets  first  bunch  near  the  grid  and  are  then 
jettisoned  at  approxima tely  equal  intervals. 

(E)  Position  of  potential  minumum  vets  us  time 

Another  presentation  of  the  system  is  that  of  the 
position  of  the  potential  minimum  against  time.  However  this 
presentation  suffers  from  an  inherent  weakness  in  the  form 
of  a  built-in  discontinuity,  since  in  our  model  there  is 
no  such  thing  as  a  true  minimum,  but  only  the  position 
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of  a  sheet  temporarily  at  the  lowest  potential.  Thus 
at  odd  moments  there  is  a  sudden  and  discontinuous 
jump  from  one  sheet  to  the  other.  Figs.  7  a  and  b 
clearly  show  the  jerky  movement  of  the  potential 
'minimum1  with  time.  In  Fig.  7  a  the  movement  is  so 

fast  thr t  we  only  get  a  blur  in  the  form  of  a  thick 
horizontal  line,  but  in  Fig.  7  b  the  jumps  are  quite 
noticeable,  especially  on  the  downward  strokes. 

(F)  Field  across  the  Ion  I  versus  time 

In  the  presence  of  the  load  it  is  often  convenient 

to  plot  the  field  E’  across  the  load  against  time,  as 

shown  in  Fig.  d.  Since  this  is  possible  only  when  the 

load  is  present,  ve  have  chosen  C  ■  9.0,  -  0.1  ns 

an  example.  This  should  be  compared  with  Fig.  6  b  where 

V  ,  is  plotted  against  t<ma  for  a  short-circuited 
mi  n 

diode  and  t  ijalso  equal  to  9»0.  In  the  case  of  Fig.  d 
the  field  E1  does  not  change  sign,  once  the  transient 
has  died  out. 

Since  our  system  is  10  in  character,  the  (x't  x') 
phase  space  is  a  20  plane  and  the  phase  trajectories 
of  the  system  are  3D  curves  with  the  t'-axis  at  right 
angles  to  the  ( x ' , x '  )-plane  .  If  we  then  'strobe'  the 
trajectories,  i.e.  note  the  values  of  x',  x'  at  regular 
time  intervals  and  plot  them  in  the  (x ', x 1 ) -plane ,  we 

obtain  what  is  usually  referred  to  as  a  Poincar£ 
section  of  the  system.  Let  us  consider  such  a  presentation 


I 
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but  at  the  same  time  alter  it  somewhat  by  taking 


x 1  ,  .  x*  instead  of  x',  x 1  . 

min  min 

(G)  Velocity  of  x*  ^  vers  is  its  position 

The  velocity  of  the  potential  minimum  x‘ 

against  its  position  is  shown  for  two  cases  respectively 

in  Figs.  9  a  and  b.  In  Fig.  9  a  we  have  L  *  4.0  and 

a  short-circuited  diode  and  in  Fig.  9  b  we  have  i  =  9.0 

and  a  load  in  the  form  of  T  »  0.1,  T  »  1.0.  N’aturally 

H  L 

this  is  another  presentation  of  the  resilts  respectively 

shown  in  Figs.  7  a  and  b.  A s  is  usual  with  other  Poincare 

sections  the  presentation  of  the  results  is  more  revealing 

and  possibly  easier  to  interpret,  since  we  now  have 

information  about  both  x*  and  x1  .In  Fig.  9  a 

.min  min 

the  velocity  of  the  potential  minimum  clearly  settles 

to  a  constant  value  i 1  .  -  x  .  /v  »  0.83  (note  that 

min  min  o 

in  fact  this  is  the  velocity  of  a  given  sheet),  its 

position  oscillating  in  a  narrow  range  near  the  middle, 

x1  *  x  .  /d  sg  0.5.  In  the  case  of  Fig.  9  b  the 

Poincar^  section  reveals  marked  variations  both  in 

the  position  x*  .  and  in  the  velocity  x'  .  . 

Regrettably  both  presentations  suffer  from  the 

artificial  discontinuities  caused  by  jumps  from  one 

charge  sheet  to  the  other  which  are  inherent  in  our 

definition  of  x'.  and  x1  . 

mm  mm 

(H)  Potential  minimum  versus  the  field 

One  more  presentation  of  the  results  proved  to  be 
of  interest.  Instead  of  plotting  against  time, 

it  is  probably  better  in  the  presence  of  a  load  to  plot 
Vmin  against  R*.  This  type  of  presentation  is  quite 


25 


smooth,  i.e.  it  does  not  suffer  from  artificial 
discontinuities,  as  can  be  seen  from  ? ig.  10  which 
is  plotted  for  .  9.0,  TR  -  0.1,  TL  -  1.0.  Fig.  io 
should  be  compared  with  Figs.  ?  b  and  9  b  which  are 
Plotted  for  similar  values  of  the  relevant  parameters. 
Both  the  periodicity  of  the  system  and  the  relaxation 
character  of  the  oscillations  indicated  by  lack  of  mirror 
symmetry  are  clearly  shown  in  the  drawing. 


The  effect  of  the  load 


'Che  majot  objective  of  these  investigations  is  to 
determine  the  effect  of  the  load  on  the  behaviour 
of  a  plasma-filled  diode,  in  particular  its  possible 
role  in  the  control  of  chaos.  So  far  we  have  developed 
a  technique  of  investigations  by  creating  a  suitable 
mathematical  model  of  the  system  and  evovlving  various 
forms  of  presentation  of  the  results.  We  now  wish 
to  describe  the  effect  of  the  load  on  the  stability  of 
operation  of  the  system  in  th*  absence  of  positive  ions 
and  then  make  suggestions  concerning  possible  influence 
of  positive  charges. 

The  mathematical  model  described  in  Section  2 
includes  an  external  load  which  is  in  the  form  of 
a  resistor  and  an  inductor  in  series  or  a  resistor 
and  a  capacitor  in  parallel.  Fig.  1,  the  irreducible 
capacitance  of  the  tube  being  connected  in  parallel 
with  the  load,  the  source  of  power  being  in  the  form 
of  a  constant  current  t  injected  through  the  grid. 

In  our  case  the  simplest  load  takes  the  form 
of  a  pure  resistor  R,  as  has  been  briefly  discussed 
elsewhere  * .  In  Fig.  11  a  we  have  a  table  of' results 
summarizing  our  investigations.  We  find  that  for  low 
values  of  TD  there  are  no  oscillations  until  the 
current  o  ^  3.0.  This  is  consistent  with  the  results 
obtained  analytically  when  the  corpuscular  nature 


of  the  electron  stream  has  been  Ignored  (hydrodynamic 

g 

approximation)  ;  we  then  find  from  curve  a  in  Fig.  7 

of  ref.  8  that  in  a  short-circuited  diode  (T_  ■  0 ,  4 '  m  1.0) 

K  o 

the  virtual  cathode  appears  also  when  •»-  ^  8.0  ,  i.e.  for 
the  same  values  of  the  current  as  those  leading  to  oscillations 
in  the  numerical  case.  This  fact  hnd  already  been  noted 
by  others  However  this  simple  relationship  is  no  longer 
valid  for  larger  values  of  .  For  example  at  T^-  8.0 

and  t  -  4.0  the  system  becomes  unstable  and  begins 
to  oscillate,  bringing  the  anode  voltage  down  to  a  mean 
value  of  4*  m  0.8,  although  in  ref.  8  for  -  4.0  the 
virtual  cathode  does  not  set  in  until  -  0.35.  Clearly 

even  a  simple  resistive  load,  which  is  absent  in  ref.  8, 
has  a  profound  effect  on  the  behaviour  of  the  system, 
probably  by  intriducing  an  element  of  feedback  in  addition 
to  that  already  present  in  the  system  and  associated  with 
the  space  charge  acting  as  a  medium  facilitating  communication 
between  the  two  electrodes  (see  also  ref.  7).  >'s  we  shall 
see  shortly  the  effect  of  the  load  becomes  more  wide  ranging 
when  reactive  components  are  included.  Figs.  11  b  -  d  show 
typical  oscillations  for  three  different  values  of  TD  , 
namely  -  0.01,  0.1  and  0.8  and  for  a  single  value  of 
the  injection  current  t  -  9»0»  The  oscillations  are  rather 
irregular,  their  amplitude  and  shape  varying  in  a  random 
fashion;  also,  as  mentioned  before,  they  have  the  character 
of  relaxation  rather  than  harmonic  oscillations,  their 
rise  associated  with  sheet  dumping  at  the  grid  (see  Figs. 11  c 
and  d)  being  much  faster  than  their  downward  stroke,  which 
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is  associated  *ith  a  slower  build- j;>  of  the  interelectrode 

space  charge.  This  is  clearly  indicated  by  the  density  of 

dots,  the  system  being  'strobed'  at  eq  ial  time  intervals. 

It  should  be  noted  ttu  t  the  amplitude  of  oscillations 

increases  by  a  factor  of  j5  as  T„  /aries  between  0.01 

H 

and  0.S;  so  does  the  frequency  of  oscillations,  but  only 

by  some  1CK.  The  corresponding  frequency  spectra  of  the 

oscillations  are  shown  in  Figs.  11  e  -  g.  Ye  find  that, 

because  of  their  relaxation  character,  the  oscillations 

are  rich  in  harmonics;  also  they  become  progressively 

'cleaner'  as  T_  increases,  their  chaotic  or  nore 
H 

accurately  random  character  being  more  pronounced  at 

lower  values  of  T  .  It  should  be  added  tuat  the  two 

R 

figures  in  brackets  respectively  refer  to  the  DC  level 

and  the  amplitude  of  the  fundamental  component,  the 

amplitude  being  normalized  ti  unity  in  the  drawings. 

Let  us  now  consider  a  load  which  has  a  reactive 

component  in  the  form  of  an  inductor  L  in  series  with 

a  resistor  R,  as  shown  in  Fig.  1,  the  resilts  of  the 

computations  now  being  collected  in  a  table  shown  in 

Fig.  12  a.  Ye  find  that  for  i-  *  4.0,  i.e.  in  the  stable 

range  of  operation  (see  Fig.  i ) .  the  addition  of  an 

inductor  to  a  purely  resistive  load  of  T_  ■  0.1  does 

R 

not  caase  instabilities.  The  situation  is  qiite  different 
however  for  larger  values  of  the  injection  current,  such 
as  L  *  9.0.  Now,  as  we  know,  the  oscillations  are  likely 
to  be  present  even  in  a  short-circ  lited  diode  and  the 
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addition  of  a  purely  resistive  load  increases  their 
amplitude,  as  shown  in  Figs.  11  b  .  d.  The  addition  of 
an  inductor  further  emphasizes  the  effect,  provided 
Tl  is  sufficiently  large*  this  can  be  clearly  seen 
by  comparing  Figs.  11  c  and  12  d,  the  effect  of  T 

R 

probably  being  stronger.  Figs.  11  d  and  12  d.  A  weak 

inductor,  Tj_  «  0.001,  has  a  negligible  effect  on  the 

amplitude  and  shape  of  the  wave,  as  can  be  seen  from 

^i@s .  11  c  and  12  o  (note  different  scales).  However 

*  0.1  and  1.0,  Figs.  12  c  and  d,  there  is 

a  considerable  growth  in  the  amplitude  of  oscillations, 

finally  reaching  a  factor  of  ~  J.  fit  the  same  time  the 

oscillations  again  become  'cleaner1,  i.e.  more  regular, 

although  they  still  retain  their  relaxation  character. 

This  is  particularly  noticeable  in  frequency  spectra. 

Figs.  12  e  -  g,  which  clearly  show  the  'cleaning  up' 

process;  at  the  same  time  the  relaxation  ch;  racter  of 

the  oscillations  is  further  confirmed  by  the  presence 

of  harmonics,  Fig.  12  g,  apart  from  which  the  spectrum 

is  remarkably  'clean1.  There  is  some  suspiscion  here 

that  we  may  be  dealing  with  a  resonance  phenomenon, 

bearing  in  mind  the  presence  of  R,  L  of  the  load  and 

the  inherent  capacitance  Cy  of  the  tube,  in  our  opinion 

this  point  merits  further  investigation,  especially  since 

some  preliminary  computations  revealed  very  high  amplitudes 

of  'clean'  oscillations  for  vanishingly  small  T  and 

R 

moderate  values  of  T^. 


and 
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Tn  our  third  aiirl  lrst  example  the  load  is  in 

the  form  of  a  resistor  R  end  a  capacitor  C  in 

parallel,  as  shown  in  Fig.  1.  In  fact  this  is  equivalent 

to  a  purely  resistive  load  with  an  abnormally  large 

capacitance  assigned  to  the  tube  and  given  by  Cy  ♦  C 

instead  of  just  C  .  However,  as  can  be  seen  from 

(2.31),  (2.31)  and  (2.34),  (2.25)  this  analogy  is  not 

complete,  since  T_  does  not  always  appear  in  the 

K 

combination  (T^  +  Tc )  -  R(Cy  ♦  C)/tQ  .  The  results  of 

our  computations  are  collected  in  a  table  and  shown 

in  Fog.  13  a .  Ye  now  find  that  in  this  case  the 

oscillations  only  occur  for  u  «  9.0,  Figs.  13  b  -  d. 

The  three  sets  of  results,  *  0.1  and  three  values  of 

Tc  ■  0.001,  0.1  and  1.0,  are  very  simil  r  and  they 

hardly  differ  from  those  for  T„  -  0.1,  T,  -  T_  »  0 

shown  in  Fig.  11  c.  This  further  supports  our  conjecture 

that  the  load  capacitance  C  is  merely  added  to  that 

of  the  tube  Cy ,  no  resonance  phenomena  being  possible 

in  this  configuration.  Comparing  Figs.  11  f  and  12  f 

we  find  that  there  is  no  noticeable  change  in  the  frequency 

of  oscillations  due  to  the  addition  of  C  ,  although 

the  frequency  spectrum  is  slightly  ‘cleaner1  in  the  presence 

of  C.  This  is  even  more  noticeable  when  we  compare 

Figs.  11  f  and  13  g  where,  in  the  latter  case,  T_  ■  1.0. 

c 

One  would  assume  that  the  R,C  combination  acts  here 
as  a  filter. 


Finally  it  is  worth  noting  the  effect  of  various 

types  of  load  on  the  intermediate,  or  hysteresis  range 

of  values  of  the  injection  current  L  ,  Fig.  J  .  ie  have 

already  mentioned  the  fret  that  in  the  case  of  a  purely 

resistive  load  oscillations  may  occur  even  for  injection 

currents  ns  low  ns  v  -  4.0,  provided  that  is 

sufficiently  large,  in  our  case  ■  0.8,  Fig.  11  a. 

The  same  applies  to  the  intermediate  range  of  injection 

currents,  say  t  -  6.0,  where  hysteresis  occurs  in  the 

case  of  r  short-ci rc  li ted  diode.  Fig.  J.  In  the  presence 

of  a  relatively  large  inductance,  T_  *  0.1,  T.  »  1.0, 

K  L 

we  find  from  Fig.  12  a  that  although  no  oscillrtions 

occur  for  b  *  4.0,  they  are  excited  at  w  »  6.0; 

however  this  only  happens  after  a  realtively  lengthy 

transient,  ns  If  the  system  could  not  quite  decide 

what  to  do.  Fig.  14.  At  the  same  time  the  hysteresis 

effect  clearly  shown  in  Fig.  J>  disappears,  the  system 

following  the  lower  branch  of  the  curve  in  both  direction 

In  the  case  of  a  capacitance  C  in  paralle  with  R  no 

oscillations  are  excited  at  (,  *  6.0,  T  -  0.1  and 

K 

Tj,  ■  1.0,  as  is  shown  in  Fig.  1?  a. 
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6.  Chaos  prolegomena 

In  the  discussion  of  chaos  it  is  always  important 
to  define  precisely  the  system  and  the  corresponding 
phase  space  and  its  Poincar£  sections.  In  our  case 
the  electrodes,  the  space-charge  which  in  general 
would  comprise  negative  electrons  and  positive  ions 
and  the  load  together  represent  our  system,  ^ince 
we  have  energy  losses  in  the  system,  both  at  the 
electrodes  where  the  kinetic  energy  of  the  charged 
particles  is  being  converted  into  heat  and  in  the  load, 
where  the  flow  of  current  genera  res  the  usual  heat  losses 
in  the  resistive  component  of  the  load,  the  system  is 
clearly  non-Hamiltonian  and  we  should  expect  the 
appearance  of  strange  attractors  if  and  when  the  system 
becomes  chaotic. 

At  this  stage  of  our  investigations  we  assume 
that  the  electrodes  are  infinite  planes  which  are  parallel 
to  each  other,  so  that  the  envisaged  system  is  ID  in 
character,  the  corresponding  phase  space  being  constructed 
from  the  position  and  velocity  variables  x  and  x  . 

The  difficulty  arises  In  deciding  which  x  and  which  x, 
or  possibly  some  other  more  s  litable  parameters  should 
be  used  in  the  representation  of  our  system. 


The  simplest  possible  model  of  our  system  is  that 
of  a  single  sheet  of  electrons  traversing  the  interelectrode 
space.  The  corresponding  equation  of  motion  for  a 
short-circ  lited  diode  has  been  derived  elsewhere  1  and 
in  the  notation  of  Section  2  is  given  by  the  following 
expression 


d*x  1 
dt  •* 


) 


'0'~(x'-f ) 


which  is  similar  to  (2.10)  for  X  *  M  ■  1  and 

2 

,ftl~  ■  2t/9.  The  solutions  of  (5.1)  are  given  by 


(5.1) 


x' 


dx  * 
dt ' 


i(l  -  cosh  m  *  t 1  )  ♦  sinh  '•) 1 1  1 

-  1  sinh  <o*tr  ♦  cosh  m't1 


(5.2) 

(6;?) 


A  set  of  the  corresponding  phase  trajectories,  each  for 
a  different  value  of  the  parameter  m '  -  21-/9  and  thus 
related  to  the  injection  current  density  c  ,  are  shown 
in  Fig.  14  a.  'Ye  find  that  there  are  only  two  possibilities: 
either  the  sheet  reaches  the  anode  (t>'<  2)  or  it  returns 
to  the  cathode  (m'>  2),  a  singular  case  being  given 
by  'o'  -  2,  it  takes  an  infinite  time  for  the  sheet  to 
reach  the  centre  of  the  interelectrode  space.  If  we  now 
strobe  the  phase  trajectory  at  equal  time  intervals  At', 
we  obtain  what  could  be  called  somewhat  grandeloquently 
a  'Poincare  section1  of  the  system.  Our  problem  in  this 
simple  ease  is  that  the  system  is  aperiodic  and  therefore 


there  is  no  preferred  value  r>f  At1,  so  that  the  whole 
point  of  plotting  a  Poincare  section  is  virtually  lost. 

Let  us  now  consider  a  more  realistic  case  when  we 

have  some  twenty  sheets  injected,  M  of  them  remaining  in 

the  interelectrode  space  at  any  given  time.  The  corresponding 

'Poincare  sections'  for  w  -  4.0  and  zero  load  are  shown 

in  Figs.  14  b  and  c,  respectively  with  and  without  the 

transient,  the  system  being  strobed  at  At'  ■  0.01,  which 

is  the  same  as  the  integration  interval.  Ve  can  see  from 

Fig.  l4  c  that  the  phase  space  trajectory  is  similar 

to  any  of  the  trajectories  shown  in  Fig.  14  a  for  m 1 < 2 , 

the  system  remaining  stable  throughout.  The  effect  of 

the  transient  is  clearly  shown  in  Fig.  14  b  -  it  simply 

amounts  to  an  upward  swing  before  the  trajectories  settle 

to  a  value  dx'/dt'  -  1  at  the  anode,  as  they  ought  to. 

The  situation  is  radically  different  for  c  -  9.0,  T_  ■  0.1, 

R 

T^  -  1.0,  vhen  there  are  well  established  oscillations. 

Fig.  12  d.  The  corresponding  phase  space  trajectories, 
again  strobed  at  At'  ■  0.01,  are  shown  in  Figs.  14  d  and  e, 
respectively  with  and  without  the  transient.  Ve  can  clearly 
see  from  Fig.  14  e  that  some  sheets  land  on  the  anode  and 
other  sheets  return  to  the  cathode,  frequently  after  some 
considerable  hesitation  in  the  vicinity  of  x'  »  l/5. 

This  is  very  much  in  agreement  with  the  time  trajectories 
shown  in  Fig.  5  b  but,  as  we  well  know,  in  total  disagreement 
with  the  simple  hydrodynamic  model  .  Clearly  the  'either, or' 
situation  with  a  singular  solution  for  m'  -  2  and  shown 
in  Fig.  14  a  also  no  longer  applies. 
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Vhen  the  system  becomes  unstable  and  begins  to 

oscillate  it  is  often  more  revealing  to  use  a  phase 

space  described  in  Section  4h  and  comprising  the  E  field 

across  the  load  and  the  value  of  the  potential  minimum  V  .  . 

min 

In  Fig.  15  a  are  shown  the  corresponding  phase  trajectories 
for  i*  ■  9.0,  Tp  •  0.1  and  ■  1.0,  which  are  the  same 
values  of  parameters  as  those  chosen  for  the  systems 
shown  in  Figs.  10,  12  d  and  14  d,  e.  In  this  case  we 
have  a  relatively  'clean*  set  of  oscillations  and  the 
plots  of  Figs.  10  and  15  a,  the  latter  including  a  transient, 
confirm  that.  Since  we  now  have  a  well  defined  period 
of  oscillations  given  by  At*  ■  2.3,  a  more  representative 
PoincarA  section  in  shown  in  Figs.  15  b  and  c,  respectively 
for  two  different  phase  values,  the  system  now  being 
strobed  once  during  each  cycle. 

Finally  in  Figs.  15  d  -  f  are  plotted  phase  space 
trajectories  for  C  »  9.0,  •  0.1  where,  as  can  be 

seen  from  Fig.  11  c,  the  oscillations  are  particularly 
'dirty*.  Ve  now  find  that  the  corresponding  phase  space 
trajectories  exhibit  similar  characteristics.  Fig.  15  d 
being  strobed  at  &t*  »  0.01  (integration  interval)  and 
Figs.  15  e  and  f  «t  it*  ■  2.1,  which  is  the  approximate 
period  of  oscillations.  Although  this  is  still  not  a 
behaviour  which  could  be  described  as  being  'chaotic', 
it  certainly  is  much  more  irregular  than  that  shown  in 
Figs.  15  a  -  c. 


In  our  searcu  lor  a  convincing  example  of  chaos  in 
the  absence  of  positive  ions  ve  have  investigated  the 
effect  of  very  high  injection  currents  in  the  case  of 
a  short-circuited  diode.  In  Figs.  16  a  -  e  ve  have 
plotted  \ 1 .  against  **  respectively  for 
L  »9,  12,  15,  18 .  anu  o.  t'e  note  that  initially  during 
the  transient  the  pattern  is  full  of  discontinuities 
or  jumps  and  then  it  settles  down  to  a  rather  fussy 
eye-shaped  curve  reminiscent  of  Fig.  9  b  and  representing 
relaxation  oscillations.  As  the  valui  >f  the  injection 
current  t,  inci  eases  >e  find  thr  t  first  a  gap  appears 
in  the  closed  curve  indicating  that  some  combinations 
of  xj^  ,  x4in  are  n0  lon£er  possible.  This  gap  tends 
to  grow,  the  fussy  eye-shaped  pattern  shrinking  more 
and  more,  until  at  t_  «  18  and  2 3  we  are  left  with 

a  fairly  shapeless  distribution  of  points  in  what,  in 
en'ect,  is  a  Poincare  section.  At  the  same  time  the 
transient  also  becomes  very  scattered.  Although  the 
pattern  has  become  very  irregular,  before  we  can  call 
it  chaotic  we  still  have  to  find  a  hierarchy  of 
bifurcations.  Ve  know  from  the  literature  '  that  in 
a  short-circuited  diode  there  is  at  least  one  'bifurcation 
point  at  C  •  8;  for  there  are  two  solutions  of  the 

corresponding  hydrodynamic  equation,  i.e.  (6.1)  with  ot  «  o 
one  stable  and  one  unstable  and  that  at  l  >  3  the 

two  solutions  coalesce.  For  i  ^  3  there  is  a  sudden 
discontinuity  in  the  current  reaching  the  anode.  Fig.  3. 
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if  we  follow  the  hydrodynamic  model  and  in  the  numerical 
calculations  the  system  breaks  out  into  oscillations. 
'Yhether  this  finally  leads  to  a  )ro,ier  chaotic  behaviour 
as  v  increases  even  further  requires  additional 
investogat Lons . 
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7,  The  influence  of  positive  ions 

So  far  we  h.  ve  only  considered  the  c  >nditions 

on  a  vacuum  diode,  i.e.  when  the  space  charge  is  solely 

due  to  the  negative  charge  of  the  electrons.  However 

our  main  task  is  to  consider  the  conditions  in  a  plasma 

filled  diode,  i.e.  in  a  system  where  the  space  charge  is 

due  to  the  negatice  ch<  rge  of  the  electrons  and  the 

corresponding  positive  charge  of  ions,  the  initial  ratio 

of  the  two  being  left  to  our  choise.  The  idea  of  deliberately 

adding  positive  ions  was  probably  first  suggested  by 

J,R, Pierce  ^  and  at  the  time  it  was  simply  intended 

to  increase  the  perveance  of  the  system,  so  that  more 

current  could  be  drawn  from  the  cathode.  Ve  now  know 

that  the  sit  iftion  is  much  more  complicated,  the  system 

being  prone  to  all  sorts  of  instabilities  and  oscillations, 

7 

including  chaotic  behaviour 

The  simplest  way  of  taking  into  account  the  effect 
of  positive  ions  is  to  assume  that  they  form  an  immobile 
backdrop  against  which  the  movement  of  electrons  takes 
place.  The  assumption  of  immobility  is  based  on  the  fact 
that  respective  masses  of  ions  and  electrons  usually 
differ  by  at  least  three  orders  of  magnitude.  If  we  then 
assume  that  both  electrons  and  ions  form  perfectly  smooth 
space-charge  clouds,  the  corresponding  differential 
equation 


J  <  0 
o 


(7.1) 
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can  be  solved  analytically  Here  Jq  is  the  injection 

'V 

current  density  rnd  4  -  m  v  "/2e  expresses  the  initial 

o  ©  o 

electron  velocity  vq  of  the  injected  electrons.  Since 

^  . 

•  -  n,  /n  expresses  the  contribution  due  to  positive 

XO  eo 

ions,  their  absence  is  indicated  by  ot  ■  0;  (6.1)  then 

3 

reverts  to  its  familiar  form  valid  for  a  vacuum  diode 
The  main  difference  between  vacuum  nnd  plasma-filled 

diodes  is  that,  as  we  can  see  from  (6.1),  in  the  absence 

2  ,  ,  2 

of  positive  ions  d  p/dx  must  remain  positive,  whereas 
in  their  presence  the  second  derivative  may  change  sign, 
for  example  when  for  sufficiently  large  4  the 

electron  density  drops  below  that  of  the  ions. 

In  the  computational  model  of  our  system  the  electron 
cloud  is  expressed  by  a  series  of  infinitely  thin  2D 
sheets  of  charge;  for  a  single  sheet  the  presence  of 
immobile  positive  ions  is  then  represented  by  a  correcting 
factor  f  ■  1  -  &  in  (2.2),  so  that  for  a  ■  1  a  single 
sheet  would  experience  no  forces  and  would  travel  with 
a  constant  velocity  from  one  electrode  to  the  other. 
Mathematically  this  can  be  expressed  by  writing  in  place 
of  (2.2)  1  : 


-tr( 


(7.2) 


In  the  presence  of  several  sheets  the  situation  is  more 
complex  since  in  place  of  (2.3)  we  now  have: 
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(1.3) 


E.  -  Btt(«11  sheets  to  the  left  of  j  )  ♦  E.  ♦ 

J  A  A  J 

E^(all  sheets  to  the  right  of  j) 

•  EIl(i>  J>  *  Ej  + 

*  -tHl^.-d1'  - j  +  0  -  -  0} 

bearing  in  mind  that  «  -  Cl  peQ  »  -  £X Np  /A.  In 

the  absence  of  positive  ions  <?  -  0  and  (7.3) 
reduces  to  (2.3). 

Turning  to  the  equations  of  motion  we  now  find, 
repeating  the  algebra  contained  in  (2.4)  -  (2.9)  that 
(2.10)  now  acquires  a  more  general  form 

-fr^-  ‘  -ftHt-l,  xl  -  j  ♦  *]  -  aN(xj  -  (  -4) 

Consequently  in  the  presence  of  immobile  positive  ions 
all  the  expressions  derived  previously  and  quoted  as 
(2.17),  (2.25),  (2.3D  and  (2.34)  are  still  valid  as 
long  as  we  substitute  everywhere  the  RHS  of  (1.4) 
for  the  corresponding  shorter  expressions  iu  the  square 
brackets.  One  can  look  at  the  new  expressions  as 
computational  models  of  (1,1)  when  a  siitable  load  is 
connected  across  the  electrodes. 


rwr”" 


At 
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In  order  to  obtain  an  expression  for  the  new 
potential  distribution  in  the  presence  of  immobile 
positive  ions  we  have  to  calculate  first  of  all  the 
additional  potential  distribution  due  to  the  ions. 
Ve  find  from  (6.?)  that  the  field  at  x  due  to  the 
ions  alone  is  f»iven  by 


C.5) 


Since  <  0  the  field  for  0  x  <  -§d  and 

0  for  ^d  < x ^ d .  Integrating  (1.5)  we  now 
obtain : 


J  *• 

^III  '  “  j  EHIdx  “  ‘^[^dT  ‘  ?X}„ 

■  a  V  x  .  ( -  l} 

2  b,  j  \  <1  t 


(7.  6) 


where  *(  *  0  at  x.  »  0  and  x.  *  d,  the  function 

TTI  J  J 

being  a  parabola  with  its  maximum  at  x.  «  d,  as  was 

J 

to  be  expected.  Adding  (7.6)  to  (2.42)  we  now  obtain: 


*i  ■  7^h<Jd-£,v  *  4^  '  ***ilxi  -  d,l  '*■>  a-: 


Dividing  both  sides  by  (2.4?)  and  adding  the  effect  of 
a  load 

*3  ’  -  -fr  *u|,  *i  *  >1  -  2E'*j  *  1  17-81 

It  should  be  noted  that  since  x!^l  by  definition,  the 

3 

additional  term  associated  with  the  positive  ions  has 
as  opposite  sign  to  the  preceding  terms  associated  with 
the  electron  cloud. 
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In  order  to  confirm  the  equivalence  of  (7.1)  and 

(7.4),  we  decided  to  consider  a  simple  case  of  a  short- 

circiited  diode  with  5  «  -  and  operating  at  two  different 

values  of  the  injection  current  t  ,  vis.  c  -5.82 

and  8.00.  As  was  to  be  expected  in  the  presence  of  positive 

ions  the  relative  values  of  the  injection  current  density 

are  much  higher.  The  results  of  our  computations  are  shown 

respectively  in  Figs.  17  a  and  b.  As  before  a  hydrodynamic 

equation  of  the  type  (7.1)  provides  us  with  both  stable  and 

unstable  solutions,  however  a  numerical  solution  of  (7*4) 

invariably  leads  to  stable  configurations  only,  provided 

we  are  outside  the  range  of  oscillatory  solutions.  Thus 

in  Fig.  25  or  rof.  7  we  have  a  plot  of  stable  and  unstable 

solutions  of  (7.1)  for  ^  -  1 ,  the  first  bifurcation  point 

occurring  at  <*/rr  *  1,  which  corresponds  to  our  /2l/?tt  -  1 

or  /T-  6.66435.  To  the  left  of  this  point  we  have  stable 

solutions  for  E  -  0,  which  corresponds  to  a  constant 
o 

potential  distribution  shown  in  Fig.  17  a,  whereas  to  the 

right  of  the  bifurcation  point  we  have  stable  solutions  for 

E  #  0,  as  is  shown  in  Fig.  17  b,  the  F.  =  0  solutions  then 
o  o 

being  unstable;  it  is  worth  noting  that  the  maximum  of  the 
potential  function  is  very  close  to  its  analytical  value 
of  1.73525.  Having  established  at  least  to  some  extent 
a  formal  agreement  between  the  results  obtained  with  the 
help  of  (7.1)  and  (7.4)  we  did  not  persue  trus  line  of 
investigations  any  further  for  reasons  which  are  discussed 
more  fully  in  the  following  paragraph. 
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In  our  opinion  the  above  numerical  model  of 

a  plasma-filled  diode  suffers  from  two  weaknesses. 

First  of  all  the  assumption  of  infinitely  heavy 

(Immobile)  positive  ions  is  unwelcome  since  in  a  way 

it  introduces  singularities  into  the  governing 

equations  of  motion.  However  what  is  possibly  even 

more  worrying  is  the  fact  thrt  although  the  electron 

cloud  has  been  discretized  by  being  converted  into 

a  number  of  charge  sheets,  the  corresponding  cloud 

of  positive  ions  has  not.  We  know  from  other 
10 

considerations  th;  t  this  is  a  dangerous  path  to 

follow.  We  are  therefore  s  iggesting  in  a  new  proposal 
that  the  work  should  he  continued  by  discretizing 
both  olottds  of  charge.  Since  the  discretization  of 
the  electron  cloud  not  only  emphasizes  the  corpuscular 
nature  of  electrons  but  also  leads  to  new  phenomena 
such  as  oscillations  which  could  not  have  been  predicted 
analytically  using  n  perfectly  smooth  model  of  the  cloud, 
it  is  to  be  expected  that  the  discretization  of  the 
positive  ion  cloud  may  well  lead  to  similar  new  phenomena 
which  could  not  be  predicted  at  present.  This  would  be 
of  particular  relevance  in  the  analysis  and  control  of 
chaos  in  a  plasma-filled  diode. 


L 
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CONCL  JSIONS 


In  the  first  period  of  our  investigations  we  have 
been  able  to  establish  s  litahle  numerical  models 
describing  the  conditions  inside  a  plane  diode  with 
and  without  positive  ions.  Having  developed  s  litable 
computer  codes  ve  have  been  able  to  investigate  in 
some  detail  the  effect  of  load  in  a  system  containing 
an  electron  cloud  represented  by  a  number  of  charge 
wheats.  Ve  have  then  derived  suitable  expressions  and 
developed  appropriate  codes  for  a  system  comprising 
electron  sheets  and  an  immobile  cloud  of  positive  ions. 

Ve  have  also  made  suggestions  for  a  further  development 
of  this  particular  aspect  of  our  investogations .  Finally 
in  Section  6  we  have  discussed  the  basics  of  chaos  as 
applied  to  our  particular  system 

Ve  have  now  reached  a  stage  in  our  investigations 
when  further  extension  of  our  work  should  be  considered. 
There  are  three  obvious  routes  we  can  follow  in  order 
to  obtain  a  clearer  and  more  complete  picture  of  the 
behaviour  of  our  system,  in  particular  bearing  chaos 
and  its  control  in  mind. 

1.  The  most  obvious  is  the  addition  of  a  resonant 
circuit  in  place  of  R,  L  in  series  or  R,  C  in 
parallel.  This  would  require  the  raising  of  the  euqations 
of  motion  from  the  second  to  the  thord  order,  a  process 
which  sould  not  create  any  special  computational 
difficulties.  Since  we  know  that  chaos  can  be  controlled 


by  Injection  of  a  small  signal  of  veil  defined  frequency, 
the  presence  of  a  load  in  the  form  of  a  reson; nt  circuit 
may  be  of  significance. 

4 

2.  Ve  are  well  placed  in  view  of  oar  earlier  work 
to  investigate  the  effect  of  initial  velocity  spread 
of  the  injected  electrons.  This  has  been  briefly 
investigated  before  *  and  we  know  that  a  velocity  spread 
can  wipe  out  oscillations  altogether  and  replace  them 
with  noise.  Again  this  type  of  investigation  would  be 

of  relevance  to  the  problem  of  chaos  control. 

3.  The  most  interesting  extension  of  our  work  would  be 
in  the  direction  of  liscretization  of  plasma.  So  far  two 
approaches  have  been  used  in  the  treatment  of  plasma. 

In  the  case  of  our  numerical  model  positive  ions  have  acted 

as  an  immovable  charged  background,  the  electron  cloud 

being  in  the  form  of  charged  sheets.  In  practice  this 

amounts  to  the  introduction  of  a  correction  factor  in  ("’.2) 

given  by  f  *  1  -  «  for  a  single  electron  sheet  and 

a  more  complex  expression  in  the  case  of  several  sheets,  l 7.4 

The  other  approach,  based  on  hydrodynamics  and  frequently 

used  in  analytical  investigations,  assumes  a  smooth 
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space-charge  of  either  sign  .  In  the  last  paragraph 
of  Section  6  we  have  indicated  in  some  detail  the  reasons 
why  in  our  opinion  both  the  electron  cloud  and  the  cloud 
of  positive  ions  should  be  discretized.  This  would  recognize 
the  corpuscular  nature  of  positive  charges,  a  step,  judging 
by  our  past  experience  \  0f  considerable  practical 
importance  in  the  modelling  of  interaction  processes. 
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CAPTIOVS 


1.  The  model  of  the  diode. 

2.  Position  and  potential  of  individual 

sheets:  (a)  t  -  4.0;  (b)  -  6.0;  (c)  u  -  9.0;  (d)t  -  9.0. 

3.  The  anode  current  ».»  as  a  function  of  the 
injection  current  l  for  a  short-circuited  diode. 

4.  Interelectrode  potential  distribution: 

(a)  u  .  6.0,  S.C.;  (b)  t  -  6.0,  TR  -  0.1. 

5.  Electron  trajectories: 

(a)  v  -  4.0,  S.C.;  (b)i«9.0,  T_  -  0.1,  T  »  l.o. 

K  L 

6.  Potential  minimum  against  time: 

(a)  t  .  4.0,  S.C.;  (b)  t,  -  9.0,  S.C. 

7.  Position  of  the  potential  minimum  against  time: 

(a)  t  -  4.0,  S.C.;  (b)  i  -  9.0,  T  „  o.l,  T  -  1.0. 

8.  The  field  K*  against  time: 

W  -  9.0,  TR  -  0.1. 

9.  Velocity  of  the  potential  minimum  against  its 

position:  (a)  (,  -  4.0,  S.C.;  (b)  c  >  9.0,  T  =  0.1,  T,  »  1.0 

K  L 

10.  Potential  minimum  against  the  field: 

^  m  9.0,  T  *  0.1,  T,  -1.0. 

A  M 

11.  System  performance  for  different  values  of  i  and  T  : 

R 

(a)  Table  of  results;  potential  minimum  against  time 
for  i  -  9.0:  (b)  Tr  -  0.01,  (c)  TR  -  0.1,  (d)  TR  -  0.8; 
frequency  spectrum  for  i  -  9.0:  (e)  T_  -  0.01, 

(f)  Tr  -  0.1,  (g)  Tr  -  0.3. 
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12.  System  performance  for  TR  -  0.1  and  different  values 
of  t  and  T^. 

(a)  Table  of  results;  potential  minimum  against  time.c  •  9.0: 

(b)  Tl  -  0.001.  (c)  Tl  -  0.1,  (d)  Tl  »  1.0; 
frequency  spectrum,  u  ■  9.0: 

(e)  Tl  -  0.001,  (f)  Tl  -  0.1,  (g)  Tl  -  1.0. 

13.  System  preformance  for  TR  -  0.1  and  different  values 
of  v  and  Tj, . 

(a)  Table  of  results;  potential  minimum  against  time,  t  ■  9.0: 
<b)  Tc  -  0.001.  (c)  Tc  -  0.1,  (d)  Tc  -  1.0; 
frequency  spectrum,  w  -  9.0: 

(e)  Tc  -  0.001,  (f)  Tc  -  0.1,  (g)  Tc  -  1.0. 

14.  Phase-space  trajectories:  x‘  against  x1: 

(a)  for  a  single  sheet,  S.C.;  for  many  sheets: 

(b)  L  -  4.0,  S.C.,  tc )  the  same,  transient  removed; 

(d)  t-  «  9.0,  TQ  *  0.1,  T.  -  1.0,  (e)  the  same,  transient 
K  1* 

removed . 

15.  Phase-space  trajectories:  potential  minimum  against  field: 

(a)  i,  »  9.0,  T  -  0.1,  T.  =  1.0,  At'  -  0.01; 

K  L 

(b)  the  same,  i  t'  «  2.3;  (c)  the  same.  At'  -  2.3,  different 
phase;  (d)  u  -  9.0,  TR  ■  0.1,  At*  -  0.01,  (e)  the  same. 

At'  -  2.1,  (f)  the  same.  At*  -  2.1,  different  phase. 

16.  Poincare  sections:  x'  against  x' 

min  ”  min 

(a )  i.  -  9.0,  ( b )  t,  -  12.0,  (c)  l  -  15.0,  (d)  c  -  18.0, 

(e  h m  23.0. 

17.  Potential  distribution,  s.c.,  0<  -  1 : 

(a)  1  -  28.33,  (b)  «.  -  64.09. 
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minimum  ne.r  the  centre,  no  oscillations 

minimum  ne;  r  the  centre  when  there  ore 
10  isci  nations;  for  1.0  oscillotrons 

develop  oncl  the  hysteresis  loo,)  ;tisn,),)ei>rs 

minimum  near  the  entrance  electrode  (ftrid); 
relaxation  osci  11  ••  t ions  (fast  rise,  slow  fall) 
are  invariably  present 
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